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^ I The production of the ^Pi charmonium state, h^ at large-p^ at the Tevatron is consid- 

G^ ' ered. The colour-octet contributions to this state are found to be dominant and give 
a reasonably large rate for the production of he- This should make it feasible to look 

Ph, for this resonance in the J/ip + n decay channel. This production rate is a prediction 

OhI of NRQCD, and the observation of the h^ at the Tevatron can, therefore, be used as a 

^ \ test of the colour-octet predictions of NRQCD. 
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The production of quarkonia has conventionally been described in terms of the 
colour-singlet model [0, H, which uses a non-relativistic approximation to describe the 
binding of the heavy-quark pair produced via parton-fusion processes into a quarko- 
nium state. The heavy-quark pair is projected onto a physical quarkonium state using a 
colour-singlet projection and an appropriate spin-projection. The failure of this model 
in explaining the data on the inclusive large-p^ J/ ip production cross-section measured 
by the CDF experiment p| has led to a major revision of the theory of quarkonium 
formation. The colour-singlet model is an approximation which works well only when 
the relative velocity v between the quarks in the quarkonium state can be completely 
neglected. It is possible to include the effects coming from higher orders in a pertur- 
bation series in v in the framework of an effective field theory called non-relativistic 



NRQCD is an effective field theory derived from the full QCD Lagrangian by ne- 
glecting all states of momenta larger than a cutoff of the order of the heavy quark mass, 
M, and accounting for this exclusion by introducing new interactions in the effective 
Lagrangian, which are local since the excluded states are relativistic. Beyond the lead- 
ing order in 1/m the effective theory is non-renormalisable. The physical quarkonium 
states can then be expanded in terms of its Fock- components in a perturbation series 
in V, and it then turns out that the QQ states appear in either colour- singlet or colour- 
octet configurations in this series. Of course, the physical state is a colour-singlet so 
that a colour-octet QQ state is connected to the physical state by the emission of one 
or more soft gluons. In spite of the non-perturbative nature of the soft gluons emitted, 
it turns out that the effective theory still gives us some useful information about the 
intermediate octet states. This is because the dominant transitions from colour-octet 
to physical colour- singlet states are via El or Ml transitions with higher multipoles 
being suppressed by powers of f . It then becomes possible to use the usual selection 
rules for these radiative transitions to keep account of the quantum numbers of the 
octet states, so that the short-distance coefficient corresponding to the octet state can 
be calculated and its transition to a physical singlet state can be specified by a non- 
perturbative matrix element. The cross-section for the production of a meson H then 
takes on the following factorised form: 



where F„'s are the short-distance coefficients and On are local 4-fermion operators, of 
naive dimension dn, describing the long-distance physics. The cutoff-dependence of 
is compensated by that of the long-distance matrix elements. 

The importance of colour-octet components was seen in the phenomenology of 
P-state charmonium production at large Pt in the Tevatron data. These processes 
do not have a consistent description in terms of colour singlet operators alone 0. 
However, even for the direct production of S-states such as the J/tp or ip', where the 
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colour singlet components give the leading contribution in v, the inclusion of sub- 
leading octet states was seen to be necessary for phenomenological reasons p|. The 
long-distance matrix elements are not calculable, and a linear combination of octet 
matrix-elements have been fixed by fitting to the Tevatron data |^. Independent tests 
of the S'-state colour octet enhancement are important and recent work shows that 
a different linear combination of the same colour octet-matrix elements that appear 
in the Tevatron analysis also appears in the analyses of photoproduction [|10] and 



hadropro duct ion experiments ||TT|. These analyses provide an important cross-check 
on the colour-octet contributions. Other production modes like J/ip production in 
low-energy — e~ collisions |12[, on the Z-peak at LEP [|1^ and in T decays 



have been considered, with the purpose of determining the magnitude of the octet 
contribution. 

In this letter, we compute the large-p^ production cross-section for the ^Pi char- 
monium state, he- The production of this resonance is interesting in its own right : 
charmonium spectroscopy |15[ predicts this state to exist at the centre-of-gravity of 



the Xci Pj) states. The E760 collaboration at the Fermilab have reported the 



first observation of this resonance but this needs further confirmation. It is, there- 
fore, interesting to determine the theoretical production rate for this resonance at the 
Tevatron and to know whether there is any chance that its existence can be confirmed. 
From our point of view the production of the ^Pi is interesting for yet another reason. 
Being a P-state, the leading colour- singlet contribution is already at O(f^), and at 
the same order we have the octet production of the ^Pi state through an intermediate 
^Sq state. As we shall see later, we can infer the non-perturbative matrix elements 
for the production of he from the matrix-elements of other states. One channel which 
may be suitable for the detection of the is its decay into a J /ip + n. Armed with 
the non-perturbative matrix elements we can predict the number of J/ip + tt events 
in the total sample of J/ip^s accumulated by the Tevatron experiments, provided we 
know the decay branching fraction for he —>■ J/ip + tc. This latter number has been 
estimated from spectroscopy We believe that with these inputs, a reasonably solid pre- 
diction of the production cross-section of the ^Pi charmonium state is possible within 
the NRQCD framework. It is worth mentioning again that we are considering the 
large-pr production of this state and therefore, we are on firm grounds in using the 
factorisation relation of the kind shown in Eq. ||. 

The subprocesses that we are interested in are the following: 

q{q)+g-'S^^^+q{q), 

q + q^'Sj^^+g. (2) 

The ^5*0*^ he is mediated by a gluon emission in a -El transition. The large-p^ 
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hadronic production cross-section is given as 
da 



dpT 

^ I dy I dxiXiGa/pixi)x2Gb/pix2)- — — v~^(^^ -*'^S + ILjc). (3) 

J J ZX\ — XxG^ dt 

In the above expression, the sum runs over all the initial partons contributing to the 
subprocesses; Ga/p and Gh/p are the distributions of partons a and h in the hadrons 
with momentum fractions Xi and X2, respectively. Energy- momentum conservation 
determines X2 to be 

xiXTe~y - 2r 

X2 = — — — , (4) 

2xi — xrey 

where r = M'^/s, with M the mass of the resonance, s the centre-of-mass energy and 
y the rapidity at which the resonance is produced. 

,,./4T¥.?^, .r^?^ (5) 



The expressions for the singlet and the octet subprocess cross-sections, da/di, are given 
in Refs. Jl^ and 0, respectively. 



The production cross-section for the ^Pi state is fully specified, once we have 
specified the colour-singlet matrix element for the ^Pi state (^Oi" Pi)'j and the value 

for the colour-octet matrix element that takes the octet -"^^o state to a (^Og''{'^So)'j- 
We note that the matrix element of the singlet ^Pi state is related to the derivative of 
the wavefunction of at the origin by 

oi^CPi)) = ^\R'm'. (6) 

The Tevatron data on Xc production fixes the colour-octet matrix element which 
specifies the transition of a ^Si octet state into a ^Pj state. We would expect from 
heavy-quark spin symmetry of the NRQCD Lagrangian H] that the matrix-element for 



^5*0^^ — >• he should be of the same order as that for ■^Sf' ^ ^Pi. This is because the 
essential difference between these transitions comes through the magnetic quantum 
number so that the corrections to this equality will be of O(f^) ~ 30%. For the 
derivative of the wave-function we use a similar argument to fix it to be the same as 
for the Xc states. 

We have computed the cross- sect ions for the Tevatron energy y/s = 1.8 TeV. In 
Fig. 1, we present the ^Pi production cross-section da/dpx as a function of px, where 
the cross-section has been folded in with the branching ratio of the ^Pi state into 
J/ip + IT. We have integrated over the full rapidity interval —0.6 < y < 0.6 covered 
by the CDF experiment. For the results shown in Fig. 1, we have used the MRSD— ' 
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densities [|T9|. The parton densities are evolved to a scale Q = Mt/2. For the singlet 
matrix element, we use the value extracted from Xc decays, which is (^Oi''{^Pi)'^ = 0.32 

ig and for the octet matrix element we have (^0'/^{^So)) = 0.0098 \g. With these 
inputs, we find that the cross-section for he production (folded in with the decay 
fraction into a. J/ip and vr, which we take to be 0.5% |T^) integrated over the region 
between 5 and 20 GeV in px is quite substantial. With the 20 pb~^ total luminosity 
accumulated at the Tevatron, we expect of the order of 650 events in the J/ip + tt 
channel. Of this the contribution from the colour-singlet channel is a little more than 
40, while the octet channel gives more than 600 events. The colour-octet dominance 
is more pronounced at large-p^- The octet contribution has a flatter px dependence 
as compared to the singlet which falls off rather rapidly with pj-- The shape of the pt 
distribution is also, therefore, a testable prediction of NRQCD. 

We have studied the effect on the cross-section of the variation of the parton 
densities, the scale and the non-perturbative matrix-elements. By using GRV densities 
20| instead of MRSD-', we find that the cross-section increases by about 20-25%. 



Varying the scale from Mt/2 to 2Mt reduces the predictions by roughly a factor of 2. 
A 25% variation results if we vary the colour-octet matrix element by 30% around the 
central- value quoted above. The decay branching fraction of into a J/ip + tt could 
be as large as 1%, and if we use this instead of the 0.5% used in the above calcuations 
we could have a production cross-section which is twice as large. 

Before we conclude, a few remarks about this remarkable prediction of NRQCD 
is in order. There is now sufficient empirical evidence for the failure of the colour- 
singlet model and the need to include effects beyond this model. NRQCD provides 
a suitable way of going beyond the colour-singlet model and, as discussed, has had 
considerable success in describing a large amount of data on quarkonium decay and 
production. There however exists a different approach to quarkonium production the 
so-called colour-evaporation model or, alternatively called, the semi-local duality model 



| 2l| . The idea in this model is to relate the integral of the open-charm production 
cross-section to the sum of the resonance production cross-sections. The individual 
resonance cross-sections are not predicted but are obtained by fitting a parameter to 
experimental data. Detailed results of this model for hadropro duct ion of J / have been 
obtained recently 23 1 and the comparison with the data is as successful as that 
of NRQCD. But we must bear in mind that NRQCD is a more predictable model. In 
particular, when some of the non-perturbative parameters have been determined from 
experiment it is possible to use symmetry relations to estimate other non-perturbative 
parameters. This is true of the production process under consideration in this paper. 
We are, therefore, able to predict the number of he events that we will see at the 
Tevatron, which is not predicted in any semi-local duality based approach, where the 
undetermined parameters can only be fitted a posteriori. 

In conclusion, we find that a reasonably large rate for the production of the ^Pi 
is expected at the Tevatron, with a dominant contribution from the colour-octet pro- 
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duction channel. Heavy-quark symmetry relations allow us to infer the size of the 
non-perturbativc matrix-elements and to predict the rate for production and its 
subsequent decay into a J/ip and a vr. By looking at J/ip events associated with a 
soft pion, it should be possible to pin down the elusive ^Pi resonance at the Tevatron, 
given the large number of J/ ip events already available. We beheve that this is a firm 
prediction of the NRQCD framework and may be a useful way of distinguishing this 
from other models of quarkonium formation. 
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Figure 1: da/dpT (in nb/GeV) for he production (after folding in the decay branching 
fraction of the into a J/ip + n) at 1.8 TeV cm. energy with —0.6 < y < 0.6. The 
colour-singlet, colour-octet and the total contributions are shown. 
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